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We present a derivation of the medium dependent wave function renormalization for 
a spinor field in presence of a thermal bath. We show that, as already pointed out in 
literature, projector operators are not multiplicatively renormalized and the effect involves 
a non trivial spinor dependence, which disappears in the zero temperature covariant limit. 
The results, which differ from what already found in literature, are then applied to the 
decay of a massive scalar boson into two fermions and to the /3-decay and crossed related 
processes relevant for primordial nucleosynthesis. 
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I. INTRODUCTION 

The large amount of new and precise data on the structure of the universe, which will be provided by 
the new generation of experiments , represents a new exciting source of information for particle physics 
too. In particular, new measurements on the primordial abundance of light elements will probably provide 
a severe arena where to test new models for fundamental interactions. In this respect, the theoretical 
predictions on Big Bang Nucleosynthesis (BBN), which are considered one of the great successes of the 
hot Big Bang theory, must be refined in order to reach the same level of precision of the new experimental 
data. 

One of the necessary improvement in the BBN algorithm is to take into account that all the nucleosyn- 
thesis reactions occur in a plasma of photons, electron-positron pairs and neutrinos whose temperature 
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varies, in the relevant epoch for BBN, in the range 0.1 ~ 10 MeV (nucleons and nuclei due to their 
large mass with respect to temperature may be safely considered at T = 0). This observation led sev- 
eral authors 1^,^ to add finite temperature contributions to weak transitions like n + j/g p + and 
crossed phenomena, due to thermal radiative electromagnetic corrections. The size of these effects on 
relevant cosmological observables, like the ratio n/p, which is the key parameter for the evaluation of the 
abundances of light nuclei, as D, '^He or ^He, has been estimated to be not larger than few percent 
1^,^. However, since the new experimental data will presumably be sensible to corrections of this order 
of magnitude, the study of the thermal effects is certainly physically well-motivated. 

Although there is a general agreement in the literature on how to compute finite temperature effects 
on phase-space, vertex, mass corrections and photon emission/absorbtion ||^- |]l2| , the finite temperature 
wave function renormalization still remains an open problem. Since Refs |^,^ , where unfortunately the 
problem was not addressed in the proper way, several approaches have been proposed in literature 
|l3| . They all agree on the idea of using finite temperature Dirac spinors to obtain the corresponding 
effective projection operator. The final results, however, striking differ jl^ . 

In this paper we approach this problem in a different but straightforward way. Our results, obtained 
in section II, differ from what has been obtained in for an additional term which is due to a peculiar 
property of the heat bath. It has a different spinorial structure and spoils the multiplicative character 
of renormalization of wave-function. The presence of such a term was first recognized in [ p^ , but we 
disagree on its explicit form. We then perform a comparative analysis of the several approaches and, in 
section III, discuss some physical implications, like the change on the decay process of a heavy scalar 
boson into lepton-antilepton pairs, H l^l^ , and on /3-decay n ^ p + + Ve and crossed related 
processes. Finally, in section IV we give our conclusions. 

II. FINITE TEMPERATURE PROJECTION OPERATORS 

The ambiguity in the wave function renormalization at finite temperature is basically due to the lack 
of Lorentz covariance. The thermal bath introduces a preferred frame, as the one in which the time-like 
four vector generalizing the non-relativistic temperature parameter takes the form (/3, 0). In this rest 
frame, the equilibrium particle distribution are the usual Bose-Einstein or Fermi-Dirac expressions, with 
/3 the inverse temperature. The very definition of a particle state, which corresponds in zero temperature 
quantum field theory to an irreducible representation of the Poincare group, becomes tricky and different 
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points of view may naturally lead to quite different results for the medium dependent wave-function 
renormalization, basically due to different choices for the renormalized fields. As noticed in [p^, however, 
the non-local, momentum dependent effects of the medium imply that such a field, as a local field, could 
even be not defined. 

We will take, as in [ p^ , quite a different point of view, identifying the particle states as corresponding 
to the energy poles in the propagator. In the covariant case, for a spinor field, the propagator will always 
develop poles in the Lorentz invariant or equivalently in , the effect of the interactions being the 
shift of the pole and a change in the corresponding residue. For a particle propagating in a thermal bath, 
the interactions with the surrounding plasma spoil, in general, this property, and the location of the pole 
cannot be represented by an invariant statement. The mass shift would in fact acquire a momentum (and 
frame) dependence, which is to say that particles propagating with different speeds in the bath rest frame 
will acquire, in general, different inertia due to the interactions with the plasma (spatial dispersion). 

Notice that, unlike the zero temperature case, the effect of the wave function renormalization, being 
momentum dependent as well, cannot be absorbed via the introduction of a local countertem in the bare 
lagrangian density. It therefore represents a genuine physical effect of the medium on the propagating 
particles. 

It is not the aim of this paper to investigate the fundamental aspects sketched above, but rather to 
discuss a prescription for evaluating the wave-function renormalization effects. They cannot, nevertheless, 
be hidden under the carpet, and they render our point of view, at least with our present understanding, 
only a reasonable way of dealing with this problem. 

As already mentioned, our starting point will be to identify the particle states as corresponding to 
the energy poles in the field propagator G{E,p), which can be perturbatively evaluated taking into 
account the effects of interactions with the surrounding medium. Actually, only the poles which are 
perturbatively close to the free particle ones will be considered, and not, for example, the hole branches 
described in [p^ . The particle state will be therefore characterized by a momentum p and a new dispersion 
relation E = E{p) given by the new position of the pole. The wave function renormalization, which is 
in general momentum dependent, can then be read off by evaluating the residue of G at the pole. This 
prescription has been already considered by Sawyer ff^ , but obtaining quite different results from ours. 
We will comment on this later. The relevant quantity which appears in the expressions of scattering 
cross sections and decay rates are the projection operators on positive (+), or negative (— ), energy states 
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. The fact that they modify with respect to the free field case is of course due to the wave function 
renormahzation, and is the way this effect shows up in interaction processes. In what follows, wc will 
report our results directly in terms of these quantities. Hereafter all calculations are performed in the 
medium rest frame. 

In order to develop a general approach to the finite temperature wave function renormahzation, let us 
first consider the simple and well-known zero temperature case, where the problem can be easily solved 
using covariancc, and looking for the pole shift in f. If E = Y^\f + So is the self-energy produced by 
radiative corrections at one-loop, the propagator for a particle of four-momentum = iE!,p) and bare 
rest mass toq takes the form 

n (J. ^, (|) - + mp + So) 



where S^ — Si(p^) with i = 0, 1 are two scalar functions. Note that in the denominator of (2.1) we have 
neglected second order terms in S;. Using standard manipulations one obtains the following expressions 
for the mass shift (5m and the wave function renormahzation factor Z2, respectively, 



= = mo Si|p2^,„2 + So|p2=„2 , (2.2) 

= l+Si|,2_2-h2mg^ +2mo^ ■ (2.3) 

These results can be equally obtained by expanding the propagator around the free positive energy pole 



dS 



E ^ uop = \J\p\'^ + rriQ (for the negative energy pole one proceeds in the same way), and looking for the 
perturbed value uip. The residue at this pole will provide the wave function renormahzation factor and 
give again (2.3). Expanding S^ up to first order in E ~ cUp we get in factf] 



(,^-2S,p^-mg-2moSo) ^ ^ iE ^ .f)iE - ^^''^ 

where and denote the functions and their derivatives with respect to E evaluated at E — ujp, 
which is to say, in this covariant case, at = TOq- Furthermore, ujp stands for the shifted energy pole, 
perturbatively close to ujp 

c.^^-^c., + ^Ei-f ^So , (2.5) 



* Actually the Ei should be expanded around the still unknown value co^, but the difference is of higher order in 
perturbation theory. 
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which is equivalent to (2.2) 



By substituting (2.4) in (Ejl) one gets for E : 



+ mo) (1 + 2Si + + - + So 
G ~ — jTT--^ 5- + finite terms 

l + ^i + g^i + ^^o ;^+(mo + M , . . 
= ^— ^ ^ +fimte terms , (2.6) 



where now ^ is understood to be ^ = '^plvi ~ P'l and 5m — tooSi + So, which is again (2.2). From 
this expression we immediately read off the residue at the shifted positive pole which represents the 
renormalized projector operator on positive energy states 

A+ I ^ , "^0^/ I '^oe/^ {i> + mQ + 5m) _ (/^ + toq + Jm) 

Aj^- i + SiH Sj H So = Z2 — p . (2.7) 

Thus, as well known, in the covariant zero temperature case, preserves the same spinor structure of 
the free field projector and, as expected, the value of Z2 obtained within this approach coincides with 



the expression (2.3) 



The renormalized projector operator for negative energy states (antiparticles), Ap, is obtained from 



(2.7) by simply replacing (ojp, oj^ ) {^^pi ~ ^p) and p ^ ~p (more simply p^^ ^P^)- This is true 



since the Lorentz invariance of S^ guarantees a quadratic dependence of these quantities on E. Thus, 



Si and S^ of (2.7), evaluated at the positive energy pole, are simply connected with the same quantities 
evaluated at the negative energy pole, which occurs in the expression of A^. As we will see later, the 
situation can be quite different in presence of a medium. 

We now consider the finite temperature case. As before, we first expand the propagator around the 
shifted positive (negative) energy pole. The residue of G at the pole will automatically give the perturbed 
positive (negative) energy projector operator. In the medium rest frame, the radiative correction S of 



Eq. (2.1) takes the general form 

S = A{E, \p\) i?7" - B{E, IpI) p.7 - G{E, |p1) , (2.8) 

where A, B and C are scalar functions. Note that the covariant limit corresponds to take A = B and 
all dependence of A, B and G on momentum and energy via the scalar quantity . Using the above 
expression for S we get for the propagator G 

^ {^ + mo-AEj^ + Bp-^~G) 

~ (p2 _ 2AE^ -ml + 2B\p\^ + 2moC) ' ^ ' ' 



5 



As for Eq. ( p.4[ ), also in this case we expand A, B and C near E — cUp, and by using the same notation 
we have 

14-91 A 4^ ^-'^P 7 I 4/ _ pi mo_Pi, \ 

^ , (2.10) 



with 



< = + , (2.11) 

Wp Wp 

representing the dispersion relation for the finite temperature physical state and the new positive energy 
pole. This result can be also written as 

w« = ^\p\^ + (mo + Sm)^ ^LOp+ "^Sm , (2.12) 



p 



where, from ( 2.11 ), the mass shift Sm is given by 



2 I -1 2 

Sm = ^A-^-^B^d . (2.13) 
Too Too 



Thus for E ~ uj?- the propagator reads 



1 + 2A + u!„A' - B'^S^ ~ 



G ~ — + Too - AE-i° + i?p-7 - + finite terms , (2.14) 

with |} — ^pjQ — P'l- Using the expression of 5m and rearranging the terms we get for the residue at 
uj^, i.e. for the renormalizcd projector on positive energy states 



Kl=[l + A + ujpA' -\^B' -^C'\^*±^ + ^ 
^ ^ ^ ujp ujp J 2uj^ 2ujp 



\P\ 



P-1 

Too 



(2.15) 



with TOfl — mo + Sm. 

From this expression we read the two effects produced by the interactions with the medium. First of 
all the zero temperature mass is shifted to the new effective value mji, which is in general momentum 
dependent. The second one is a genuine contribution due to the wave function renormalization, which 
consists not only of a non trivial multiplicative factor but also introduces an additional term proportional 
to _B — A and with a different spinorial structure. This is the peculiar signature of breaking of Poincare 
invariance which introduces a spinor dependence in wave function renormalization. Notice that both 
terms are momentum (and frame) dependent. It is easily seen that at zero temperature, i.e. in the 
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covariant limit, A ^ B, the second term vanishes and the first one reproduces the usual multiplicative 
renormalization factor. 

In order to obtain the renormalized projector on negative energy states one should perform an analogous 
computation but expanding now all the involved quantities around the negative energy pole —ujp. In the 
general case, A, B and C do not have definite properties of parity with respect to E and so and 
can have quite different expressions which are not simply connected one each other. This is strongly 
related with the property of charge conjugation invariance of medium and interaction. For example, for 
e~ and in thermal equilibrium with photons one only expects their chemical potentials to be opposite. 
However, if the two distributions are equally populated, and thus chemical potentials vanish, the electron- 
positron plasma is charge conjugation invariant and hence A, B and C are even functions of E. This is 
for example the situation occurring at the time of BBN. Hereafter we will make this assumption, though 
more intriguing situations could be considered as well. Under the assumption of charge conjugation 
invariance, results to be 



^ V ' oj.p ujp J 2loR 2ujp 



mo 



(2.16) 



The results of (2.15) and (2. If) are, in a sense, half way between the approaches of Ref.s |^ and jl^ . 



The simple multiplicative factor, in fact, is exactly the one obtained in [^, but in that approach the 



additional term is absent. Denoting with A^^ the projector obtained in H we get 



B- A 



2ujp 



21 



Too 



(2.17) 



Actually the approach followed there is substantially different than ours. They start introducing finite 
temperature spinors "0, chosen as the solution of the nonlinear Dirac equation 

(;i - too - S)Vi = , (2.18) 

whose corresponding creation and annihilation operators are assumed to satisfy ordinary, zero tempera- 
ture, anticommutation relations. Expanding the propagator in terms of these spinors they obtain a wave 
function renormalization factor which is only multiplicative. We think that the assumption made on the 
canonical spinor basis to be used is responsible for this feature and represents the essential difference 
with our approach. As mentioned, we do not make any hypothesis on the renormalized field, if any, to be 
used, but simply recover the particle content and the corresponding projector operators from the poles 
of G and their residues, respectively. 



The fact that the renormahzed spinors are related to the free Dirac ones via a momentum dependent 
transformation in spinor space has been first stressed in In this analysis the projector operators 

are deduced, as in our approach, looking at the residue of the propagator, and the following result is 



obtained, in term of our 



2ujp 



21 



mo 



(2.19) 



It is still unclear to us the reason for the difference between and our result, since the starting point, 
in both cases, is the same. It should be stressed, however, that does not reproduce the expected 
behaviour in the zero temperature, covariant limit. In fact, since we have already shown that A^ gives 
back in this limit the correct result, A^ will still contain an explicitly noncovariant term, only vanishing 
if C = mo-B', which is not guaranteed by any general principle. 

III. APPLICATIONS TO SIMPLE PROCESSES 

In order to show the physical differences of the three different approaches to thermal wave function 
renormalization, simply summarized by the use of projectors A^ , A^^ and A^, let us discuss a relevant 
example like the leptonic Higgs decay H l^l^ . A similar analysis has been performed in At 
tree-level the matrix element for this process reads 

A{H^in-)^-'^u{l-)v{l+) , (3.1) 

where m; stands for the zero temperature lepton mass, and u(l~), v{l~^) denote the T = Dirac spinors for 
lepton and antilepton, respectively. In a thermal bath of equally populated lepton-antilepton pairsf], the 
above decay width acquires an additional contribution, due to the lepton wave function renormalization, 
which at one loop results to be 

\M\l^ ^ \A{H ^l+{X+)l-i\-))\' ^ ^TT[A+il-)A-{l+)] 
x+.x- ^ 
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g^mf 

4M2: 

I 



m 



Tr [(A+(r) + SA+{1-)) (A-(/+) + SA-{1+))] 
{Tr [Ao+(/-)Ao (/+)] + Tr [<5A+(r )Ao-(/+) + A+ {l-)6A- (1+)] } 



|7W|J + (5|X|^ ■ (3.2) 



^We assume the Higgs mass much larger than the lepton temperature so that the decaying particle can be 
considered at T = 0. 
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In Eq. (3.2) A* stand for the lepton polarizations and with A^(ZT) and SA^{P) we denote the zero 



temperature lepton projector and its first order thermal correction, respectively. The second term in the 



r.h.s. of (3^), S \M\ , then represents the thermal contribution to the decay width due to thermal wave 



function renormalization. The expression for 5A^{P) can be extracted from of Eq.s (2T5),(2T^), 



or Ag. Once substituted in A |Al| we get in the cm. reference frame 



6 M 



R 



7 ? 



5\M\ 



2 2 



Mf 
A/2 



2A-B + —A' 1 

2 2 



Ml 



'2mi g,' 
Mh 



A+^iA'-B' 



imf\^, .... 



H 



(3.3) 
(3.4) 
(3.5) 



We have not included the effect due to the renormalization of the lepton mass, in order to disentangle 
the two different effects. To consider its contribution it is sufficient to substitute the value of to/ and 
bjp with mf = mi + 5m and oj^ = ojp + Suj reported in ( 2.13 ) and ( 2.12 ) respectively, in the tree level 
squared amplitude 



\M 



2M2, 



Ml 



(3.6) 



A similar shift can also be performed in any of the contributions (3.3)-( |3.5| ), since the difference we 
introduce in this way on the squared amplitude is of higher order in perturbation expansion. 

All approaches give the expected behaviour for the squared amplitude^, which must vanish if Mh = 



2m/. Nevertheless, in the expressions (|3.3| )-(3.5) the thermal corrections take different forms. It is worth 



observing that both our result and the one obtained using A^^ coincide in the covariant limit and give 
the expected Z2 factor reported in (pT^). This is not the case if A^ is instead used. 

Another example where to compute the additional contribution due to the finite temperature renor- 
malization wave function for electrons is the /?-decay n — > p + + and the crossed related processes. 
Actually, the need for a consistent way of computing thermal corrections to these processes, relevant for 
primordial nucleosynthesis, represents one of the main motivation for our present analysis. According to 
the notation adopted in a simple computation shows that, in the nonrelativistic limit for both neu- 
tron and proton, the use of the renormalized projector Aj for the outgoing electron gives the additional 
contribution to the rate per neutron 



*In this respect we note that the result quoted in |14| as equation (11) is incorrect. 
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2tt-^ Jq \ Up ujp J 

X [1 - F,{Q)){1 - F^iujp)) , (3.7) 

where p is electron momentum, tOp its energy, Q = Af„ — Mp — LUp the neutrino energy and F^ e are 
the neutrino and electron Fermi-Dirac distributions, respectively. The corresponding corrections due 
to electron (positron) wave-function renornialization for the other crossed processes are easily obtained 
replacing the thermal factors and the expression of Q as reported in Table 1 of rcf. j|] . Despite of the 
different forms of A^, and A^^, it is interesting to note that unlike the scalar boson decay previously 



considered, all give the same result (3.7) in the nonrelativistic limit. The reason is that the extra pieces 



which make different the three expressions (see Eq.s ( p7| ) and ( |2T9| )) give a correction which vanishes 



after the angular integration of the differential rate is performed. This fact has been already noticed in 
0- 

IV. CONCLUSIONS 

In this paper we have addressed the problem of finite temperature renormalization effects on spinor 
wave function. This issue has been previously considered by many authors, with different results. The 
main question seems to be whether this renormalization effect is simply multiplicative, as in ordinary 
quantum field theory, or rather involves a non trivial spinorial dependence as well. Our starting point has 
been to look for the residue at the energy poles of the propagator, corrected for thermal interactions with 
the surrounding medium. These poles give the new energy-momentum dispersion relation and illustrate 
the particle content of the theory. This point of view bypass the difficulty and ambiguity related with 
the choice of the renormalized fields to be used at finite temperature, which are due to the non-local 
and frame dependent character of the medium effects. The projector operators on positive and negative 
states we have obtained show that, as pointed out in the breaking of Poincare invariance leads 
to the appearance of extra contributions with a different spinorial structure. They reduce to ordinary 
renormalized projectors in the covariant limit. 

The physical interest for this problem is mainly connected with the study of thermal corrections to 
nuclear interactions in the early universe, during the primordial nucleosynthesis epoch. Notably, our 
results for the case of /3-decay and related processes completely agree, in the non relativistic regime, with 
the ones which can be obtained using different approaches or results for wave function renormalization. 
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since the extra contributions mentioned above play no role in this limit. This is not the case for the decay 
of a scalar boson in fermion-antifermion pair where we have found using our approach, a different result 
for the corresponding change in the decay rate. 
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